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1. Summary and Introduction 


The purpose of this manual is to aid in the design of an aerodynamics test of an Earth or 
planetary entry capsule in a ballistic range. In this manual, much use is made of the results and 
experience gained in 50 years of ballistic range aerodynamics testing at the NASA Ames 
Research Center, and in particular, that gained in the last 30 years, while the author was 
working at NASA Ames. The topics treated herein include: 


¢ Data to be obtained; predicted flight data needed to design test 

¢ Reynolds number and dynamic similarity of flight trajectory and ballistic range test; 
capabilities of various ballistic ranges 

¢ Calculations of swerves due to average and oscillating lift and of drag-induced velocity 
decreases 

¢ Model and sabot design; materials, weights and stresses 

¢ Sabot separation 

¢ Launches at angle of attack and “slapping” with paper to produce pitch/yaw oscillations 
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Nomenclature 


area 
acceleration 

drag coefficient 

lift coefficient 

pitching moment coefficient 


pitch damping coefficient 

derivative of pitching moment coefficient with respect to angle of attack 
yawing moment coefficient 

yaw damping coefficient 


side force coefficient 
speed of sound 
diameter 


fractional velocity drop over distance of one pitch oscillation wavelength 


Young's modulus of elasticity 
strain energy per unit volume 


strain energy per unit volume averaged over shell thickness 
strain energy per unit mass averaged over shell thickness 


strain energy per unit mass averaged over shell thickness and averaged over 


entire shell 
force 


moment of inertia for pitching motion about center of gravity for model or flight vehicle 


angular impulse 

moment of inertia of shell about neutral axis per unit length 
linear impulse 

constant, see Eq. (4) 

kinetic energy per unit mass 

lift 

length of range 

Mach number 

average Mach number 

bending moment 

mass 

paper mass loading (mass per unit area) 
paper mass loading (mass per unit area) 
period of pitch oscillation 

range pressure 

= ps/p 

Reynolds number 

radius 

radius of gyration 

temperature 

sideways swerve 


t time 

ta shell thickness 

U velocity 

U, velocity of thin shell or aerobrake with respect remainder of model 
V volume 

X distance 

Xcg_ X location of center of gravity 
y swerve to side 

Ymax Maximum swerve to side 

yY acceleration to side 

Zeg_ Zlocation of center of gravity 


Subscripts: 


cg center of gravity 

f flight 

max maximum 

r range 

r for one pitch oscillation wavelength 


Greek: 


pitch angle of attack 

rms pitch angle of attack 

yaw angle of attack 

rms yaw angle of attack 

strain 

pitch angle of attack 

maximum pitch angle of attack 
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maximum rate of change of pitch angle of attack 
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wavelength of pitch oscillation 
viscosity 

3.14159... 

gas density 
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shell density 


model or vehicle density 
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range model density 


co 
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maximum compressive or tensile stress 


a oO 


yield material yield stress 
2. Description of a ballistic range; types of testing covered in this 
manual; data to be obtained 


Figure 1 shows a sketch of the NASA Ames HFFAF (Hypervelocity Free Flight Aerodynamic 
Facility) ballistic range, which is representative of enclosed test section ranges. Key elements of 
the range are the launch gun, dump tank (to contain the muzzle blast), sabot stripper (to prevent 
the sabot fingers from flying down the range) and test section with spark shadowgraph stations 
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and cameras. The Ames HFFAF range has 16 orthogonal stations and 32 cameras. In the 
HFFAF, the shadowgraphs are shuttered by 40 nanosecond Kerr cells. Open air ranges, 
operating at atmospheric pressure, do not use a dump tank or an enclosed test section. Models 
are typically launched inside multi-piece carriers, called sabots. Figure 2 shows a model to be 
launched by a four piece sabot, comprising four sabot “fingers”, locked together in the gun 
barrel by serrations. On exiting the gun muzzle, the sabot fingers are meant to be separated 
sideways from the model by aerodynamic forces acting on the fronts of the fingers and/or 
muzzle blast forces acting on the rears of the fingers. With sufficient lateral sabot finger 
separation, the fingers will impact on the outside of the sabot stripper cone and the model will 
travel down the range unaccompanied by the fingers. 

The launch gun can be a cold gas gun for low velocities (up to ~1.3 km/s) or a powder gun 
for velocities up to ~2.5 km/s. For higher velocities a two-stage light gas gun using helium or 
hydrogen (for the highest velocities) is often used. Figure 3 shows a representative two-stage 
light gas gun. The piston (often plastic) is accelerated by the gunpowder gas to velocities 
between (typically) 500 to 1000 m/s, and greatly compresses and heats the light working gas. 
When the light gas pressure reaches a certain value, the diaphragm breaks and the 
compressed light gas accelerates the launch package (with the model) down the launch tube. 
Muzzle velocities up to 7 km/s can be obtained fairly routinely (with hydrogen). With special gun 
designs and operating conditions and light projectiles, launches at velocities up to ~11 km/s 
have been achieved. 

The vehicle models treated herein are Earth and planetary entry capsules, including those 
with low density aerobrakes (i.e., inflatables, deployables) outside a compact centerbody. In 
general, the models are axisymmetric or have a plane of symmetry. The models may trim at an 
angle of attack (and hence, lift) due to an off-axis weight or a trim tab. More complex models, 
such as airplane shapes or capsule models with an abort rocket tower attached, are not 
discussed. Also, aerodynamic heating models, which are generally launched with minimal pitch 
and yaw angles, are not treated. Testing of highly lifting and swerving models, such as those for 
shuttle tank foam debris, is also not discussed. The exclusion of a number of types of testing 
keeps the scope of the present paper within manageable bounds, and permits us to examine in 
some detail testing that does lie within the scope of the paper. 

The data sought are generally: 


The drag coefficient, Cp, as a function of the pitch angle of attack, a 

The lift coefficient, CL, as a function of the pitch angle of attack, a 

The pitching moment coefficient, Cm, as a function of the pitch angle of attack, a 
The pitch damping coefficient, Cy , as a function of the rms amplitude of the 


pitch oscillations, & 
The above data are often sought for a number of Mach numbers in a test entry 


For non-axisymmetric models with a plane of symmetry, other coefficients, such as the side 
force coefficient, Cy, the yawing moment coefficient, C,, and the yaw damping coefficient, C, 
can also be obtained and all coefficients are, in principle, functions of the pitch angle of attack, 
a, and the yaw angle of attack, B. In addition, the damping coefficients are functions of the rms 
angles of attack, @ and B. However, in general, the model behavior is not very different with 


respect to pitch or yaw oscillations and designing of a test can be done by referring only to Cop, 
C_ and Cm as functions of the pitch angle of attack, a, at the Mach numbers of interest. 


3. Data needed for test design 


In order to design a test that simulates full-scale flght, certain data are required. This 
includes vehicle physical and aerodynamic data (the latter would be estimates from CFD 
solutions) and the flight trajectory data. These will be used to calculate the degree to which 
dynamic similarity and Reynolds number similarity can be achieved between the flight vehicle 
and the range model and to estimate the model swerves due to oscillating and average lift. 
Similarity is treated in sections 5 through 7. Swerve due to oscillating lift and the effect of its 
amplitude on the accuracy of the determination of the lift coefficient as a function of the pitch 
angle of attack is discussed in section 8. Swerve due to average lift for lifting vehicles can 
obviously lead to the model flying out of the field of view of the cameras and is discussed in 
section 9. 

To perform the required calculations, considerable data, including predicted flight data, is 
needed. First, one needs the vehicle physical data, as follows: 


1 
¢ Diameter, D (area, A = zP) 


¢ Mass, m 
¢ Moment of inertia for pitch oscillations, | 


- Radius of gyration, rg, defined by 7, = I/m 


¢ Volume, V 
¢ Density of flight vehicle, ps = m/V 
¢ Location of CG, Xeg, Zeg 


¢ Vehicle shape (drawings) 
¢ Note that subscript f will denote flight vehicle parameters and subscript r will 
denote range model parameters 


Next, one needs vehicle aerodynamic data, usually computational (CFD) values at the flight 
conditions of interest. For each Mach number of interest, these comprise: 


¢ The drag coefficient, Cp 
¢ The lift coefficient, C_, as a function of the pitch angle of attack, a 
¢ The pitching moment coefficient, Cm, as a function of the pitch angle of attack, a; 


from this data, the derivative Ci: can be calculated 


Finally, one needs the trajectory data at each point along the trajectory to be tested. These 
comprise: 


e Atmospheric composition 

¢ Velocity, u 

¢ Density, p 

¢ Temperature, T 

¢ From the temperature, T, and the atmospheric composition, the speed of sound, c, 
and the viscosity, u, can be obtained 

¢ The Mach number, M, follows from M = u/c and the Reynolds number, Re, follows 
from 


Roses 
u 


4. Calculations which need to be performed 


For each point along the flight trajectory to be simulated, one must calculate Re, p,/p and A/D, 
where: 


Re = Reynolds number 

ps = flight vehicle mean density 

o = atmospheric density 

dX = wavelength of pitch oscillations 
D = flight vehicle diameter 


The three aforementioned quantities must be matched, as closely as possible, between the 
flight conditions and the ballistic range conditions in order the achieve Reynolds number 
similarity and dynamic similarity. (Further discussion is presented in Sec. 5.) Then, the 
capabilities of the ranges in question need to be plotted in the Re — p,/p plane. These 
capabilities are determined by various limits discussed at a later point. The flight trajectory must 
then be plotted on the Re — p,/p plane. This will aid one to select the most suitable range in 
which to perform the desired test campaign. At this point, it may already be necessary to accept 
some compromises, that is, to accept some mismatches in Re and p,/p. For example, it may 
prove impossible to simulate simultaneously high values of Re and low values of p/p. An option 
in this case would be to accept the Re mismatch and simulate accurately the flight values of 
p;/p in the range. 

From A/D and the spacing of the camera stations, D needs to be selected to put 6 to 8 
camera stations in each pitch oscillation cycle. Then, the gun should be selected to make this D 
value as close to 80% of the gun bore diameter as possible, without exceeding the 80% figure. 
Some compromise may be necessary here. The 80% figure is good general rule, but saboted 
models have been launched that had diameters up to 90% of the gun bore diameter. In such a 
case, there would be increased concerns about sabot failure due to the launch stresses. 

If the range is an enclosed, variable pressure range, the range pressure needs to be 
selected to match Re, if possible. If the range is an open air range, the range pressure is always 
atmospheric, which determines what Re values will be obtained. ps should be chosen to match 
(Ps/P)right. The model material then needs to be chosen to match ps as close as is possible. Multi- 
material models or models with internal hollows may be necessary to get a good match of ps. 

These calculations may need to be iterated as required. As stated previously, it may be 
necessary to accept some mismatches in Re and p,/pin the final model design and test 
condition selection. After the main design work outlined above has been done, several checks 
and additional calculations should be made, as follows. 

The normalized drop in velocity per pitch oscillation cycle and for the length of the range 


and on , where du is the velocity drop and u Is the 


r range 


average velocity. To maintain a “constant” test condition in the range (as much as possible), it 


should be checked. These are du 
Uu 


would be desirable if au were less than 0.20, but this may be difficult to achieve for a 


range 


lightweight ballistic range model encountering a relatively high density test gas. There is also a 


lower limit to a connected with having acceptable accuracy in determining the drag 


range 
coefficient. To have reasonable dynamic similarity for drag slow down between the flight vehicle 


and the model, oo should not be too different for the two cases. A difference of up to perhaps 


u 


r 
0.06 should be acceptable. 

The swerves due to oscillating a should be checked; these must be larger than a certain 
minimum value if C, , the slope of the lift coefficient with varying a, is to be determined 


accurately. 

The swerve due to average lift should be checked (only applies to lifting models). Verify that 
model will remain in view of the cameras. 

Stress calculations for model and sabot need to be performed. 


5. Calculating range simulation capabilities — example 1 


To achieve similarity between the flight vehicle and the ballistic range model, one must have 
the same Reynolds number (Re) for both and the model motion must be dynamically similar to 
that of the flight vehicle. The dynamical similarity comprises dynamical similarity of pitching (and 
yawing) motion and of the slowing down of the model and vehicle due to drag. The pitching 
motion will be dynamically similar if A/D for the model and the flight vehicle are equal, where AJ is 
the wavelength of the pitching oscillation and D is the diameter of the model or flight vehicle. 


The slowing down due to drag will be dynamically similar if a for the model and the flight 
U 
r 

vehicle are equal, where u is the flight velocity and du is the change in flight velocity over the 
wavelength 2. In general, much more importance is given to pitching motion similarity than to 
drag slow-down similarity. Also, if pitching motion similarity is achieved, drag slow down 
similarity will not be too far off if r,/D is not too different between the model and the flight vehicle, 
where rg is the radius of gyration. Thus, in general, for ballistic range work, similarity between 
the model and the flight vehicle consists of matching each of Re and A/D between the two 
cases. 

The steps to make these calculations are given below. First, calculate 4/D for flight vehicle 
for each trajectory point to be tested using Eq. (1) below. [Equation (1) was derived from Eq. 
7.140 of Ref. 1 by suppressing the minus sign of Cn, and replacing the reference dimension | by 


the model diameter, D.] 


Sar 


ome ee 
D \ pAD°C,, 


(1) 


Equation (1) gives the normalized pitch oscillation wavelength in terms of the pitching moment 
of inertia of the vehicle, |, the gas density, p, the frontal area of the vehicle, A, the diameter of 
the vehicle, D, and the derivative of the pitching moment coefficient with respect to angle of 
attack,C_ . At this point, a rough, preliminary model design is necessary. This design would 


have the same shape, as close as is feasible, as the flight vehicle. The diameter could be set to 
80% of that of the available gun, as a first cut. A material must be chosen for the majority of the 
model volume. This could be steel or aluminum or some other material. The model may need to 
be weighted to have the center of gravity in the location corresponding to that of the flight 
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vehicle. Note that both the model diameter and the model material may change during the 
2 


r 
design process. For the preliminary model design, calculate 2 Then, the relation between 


i/D and p./p for the model is calculated from 


or 
A 
0 
8? Vor 
where K = = ia (4) 


(Note that p,; for the model in Eqs. (2) and (3) is not the same as p, for the flight vehicle.) Since 
the flight vehicle and the model are assumed to be geometrically similar, V/AD is the same for 
both. Equation (3) is then rearranged as 


ae (5) 


The Reynolds number, Re, can be written as 


Re= puD _ pMcD 
u u 


By manipulating equations (5) and (6), it can be shown that 


2/3 
I | “ates 


Le 
0 R23 Reu 


Equation (7) defines two limiting lines for the facility in the p,/p — Re plane. One line is for the 
maximum value of Ap, and the other for the minimum value of Aps. For the NASA Ames 
Hypervelocity Free Flight Aerodynamic Facility (HFFAF), 2 to 2.5 cycles in the 75 ft (22.86 m) 
long range has proven optimal. Two complete cycles are needed so that at least 4 extrema of 
the pitch and yaw oscillations can be seen at the 16 orthogonal shadowgraph stations. With four 
extrema, lines can be drawn joining the 2 positive extrema and the 2 negative extrema and the 
difference between the slopes of these lines allows the damping (or growth) of the oscillations in 
pitch and yaw to be accurately determined. (An example of a curve of pitch oscillation 
amplitudes varying down the range and the lines joining the extrema is shown in Fig. 4.) In 
addition, it is desirable to have at least 6 data points per cycle, which requires a minimum A of 
9.14 m (since the spacing of the shadowgraph stations is 1.524 m). (Six data points per cycle 
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are the minimum to reasonably accurately define the amplitude and wavelength of the 
oscillation. Seven or eight points would be better, if possible.) Thus, tests in the HFFAF are 
designed so that A ranges between 9.14 and 11.43 meters. 

Typical model densities range between 0.94 g/cm? (polyethylene) to 13 g/cm’. The latter 
figure is based on machinable tungsten alloy (with a density of 17 g/cm*) with enough margin to 
locate the CG correctly. An upper limit can be defined using the maximum, reasonable values of 
dX. and ps and, likewise, a lower limit can be defined using the minimum values for these two 
variables. To exceed the upper limit would require a p; above 13 g/cm? (which is very difficult to 
achieve) or by increasing A, so that there would be less than 4 extrema of the pitch oscillations 
during the model's flight within the range test section. To go below the lower limit requires a pg, 
less than 0.94 g/cm’, or a reduced A, so that there are fewer than 6 data points per pitch 
oscillation cycle. 

We will use the supersonic flight test of the Low Density Supersonic Decelerator (L.DSD)* 
project as an example to illustrate calculating the test simulation capability of a ballistic range. 
The LDSD flight tests were to demonstrate the feasibility of a Supersonic Inflatable 
Aerodynamic Decelerator (SIAD) for Mars entry. The SIAD is an inflatable torus attached to the 
outer rim of a capsule-like atmospheric entry vehicle. Figures 5 and 6 show the vehicle and the 
ballistic range model, respectively, in the SIAD-deployed configuration. Note that the ballistic 
range model is solid metal with shape representative of the inflated SIAD. For the flight vehicle 
at Mach 3.96, the range capability limits are shown in Fig. 7 as the two slanted limit lines with 
slopes Re a (p;/p) '”. Note that the LDSD is referred to as SIAD and SIAD-R in Ref. 2. 
However, not all of the area between these two diagonal limit lines is accessible to the facility 
since there are several other limiting factors. If the smallest acceptable 4 is set at 9.14 meters 
and maximum model diameter is D = 3.56 cm (for models launched from the Ames 44.5 mm 
gun) and 4.88 cm (for models launched from the Ames 61 mm gun), these define limiting values 
of A/D. Using Eq. (5), these translate into limiting values of p./p ~ 20,000 and 10,600, 
respectively. These limiting values appear as the vertical lines defining the leftmost sides of the 
facility capability zones in Fig. 7. One could conceivably extend the limit to the left by increasing 
D. This could put the integrity of the sabot at risk, since there would be very little sabot material 
left between the model and the gun bore. Likewise one could also extend the limit to the left of 
this line by decreasing A. However, this would result in less than 6 data point per pitch 
oscillation cycle. There is another p,/p value limit on the leftmost side of the facility capability 
zone. Since the maximum allowable test section pressure is 1 atm in air (Tests in other gases 
are typically done at sub-atmospheric pressure levels.), the maximum gas density is 1.194x10° 
g/cm®. Using a minimum pg of 0.94 g/cm?® (polyethylene) yields a p/p of ~800. This limit lies well 
to the left the p,/p values obtained previously from the A/D values and thus, does not limit the 
facility capability. 

Similarly, using the largest acceptable 4 value of 11.43 meters and a minimum acceptable 
model diameter of 1.27 cm (based on acceptable resolution in the shadowgraphs). These 
define a limiting value of A/D which translates into a limiting value of p;/p of 2.45 x 10°. This 
limiting value of p./p is shown as the vertical line to the right of the facility capability zone in Fig. 
7. One could move to the right of this line by reducing D further. However, this could result in 
decreased accuracy of the data obtained, since the ratio (Shadowgraph reading error)/(model 
diameter) would be reduced. One could also move to the right of this line by increasing 2, which 
would, however, result in having less than 4 extrema of the pitch oscillations within the range. 

Another limit is the maximum possible Reynolds number (for a given Mach number) based on 
the maximum facility pressure and the maximum model diameter. For the 4.88 and 3.56 cm 
LDSD models with the inflatable deployed at Mach 3.96, these limits are shown as the 
horizontal lines near the top of the facility capability zone in Fig. 7. It is seen that the line for the 
4.88 cm model truncates slightly the facility capability zone, while the line for the 3.56 cm model 
does not affect the corresponding facility capability zone. For a one atm maximum pressure 
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range like the HFFAF, one may move above the maximum Reynolds lines shown in Fig. 7 only 
by increasing D, which can be done only by 10% or so, and even that at the risk of 
compromising sabot integrity, as discussed previously. 

At the topmost part of the facility capability zone, still another limit is encountered. From the 
definition of the Reynolds number, we write 


Re = ——+— (8) 
u (p,/ p) 
Rearranging, we get 
Pp, U 
- 5 Re 9 
p, aD (9) 


Near the top of the facility capability zone defined by the leftmost vertical line and the topmost 
slanted line discussed previously, taking D = Dmax, it is easy to show from Eq. (9) that there is a 
small region where the model density exceeds the maximum allowed density. This region is 
removed from the facility capability zone by adding the slanted lines with slopes Re «a (p,/p). 

The various limits which may bound the facility simulation capability zone in the Re — p,/p 
plane are enumerated below. Not all the limits affect the capability since some limits may fall 
outside other, more restrictive, limits. 


(1) Two lines defined by Eq. (7) with slope -1.5 in the Re — p,/p plane. One line is for 
(Aps)max and the other for (Aps)min. The facility capability lies between the two lines. 


(2) Two vertical lines determined by (A/D)max and (A/D)min. The facility capability lies 
between the two lines. 


(3) One vertical line defining the minimum p./p by pmax aNd ps min. The facility capability 
lies to the right of the line. 


(4) One horizontal line at the maximum Re defined by Dmax and the maximum facility 
pressure. The facility capability lies below the line. 


(5) One line defined by Eq. (9) with slope -1 in the Re — p,/p plane. This limit depends on 
the maximum value of ps. The facility capability lies below the line. 


The solid blue lines in Fig. 7 define the facility capability zone for Mach 3.96 and a maximum 
model diameter of 3.56 cm. For a maximum model diameter of 4.88 cm, the facility capability 
zone extends further upwards and to the left, enclosed by dashed blue lines. For the Dmax = 3.56 
cm case, the limits which appear in Fig. 7 are, from the above list, (1) (both lines), (2) (leftmost 
line) and (5). For the Dmax = 4.88 cm case, the limit (4) also appears. 


We next place the flight trajectory points on the p,/p — Re plane of Fig. 7. For each trajectory 
2 


r 
point A/D is calculated using Eq. (1). Then, for the preliminary model design, 7? is determined 


for the model. K is then calculated from Eq. (4), p;/p from Eq. (5) and the Reynolds number from 
Eq. (6). The red data points show the flight trajectory data points for the Mach number range 
2.37 — 3.96. It should be pointed out that the p,/p values for the trajectory points plotted on the 
facility capability plot of Fig. 7 are not the same as the actual flight p,/p values. This is because 
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2 Z 
§ 


r r 
D2 is different for the flight vehicle and the model. (For the flight vehicle, a is 0.0421 and for 


2 pa 
2 


r 
the preliminary model design considered herein, a is 0.0599 — the difference being due to the 


light inflatable surrounding the centerbody of the flight vehicle. This is very difficult to simulate in 
the ballistic range, since the model, shown in Fig. 6, is of solid metal construction.) The priority 
is given to having the correct dynamic similarity of the pitch oscillations (i.e., equal A/Ds) rather 
than the correct dynamic similarity of the drag slow down, which would be obtained if the 

p;/p values were matched between the flight vehicle and the model. (See discussion of velocity 
drops presented in Sec. 7.) For the flight trajectory under discussion here, and the range 
simulation of the same, the key parameters are given in Appendix A to enable the reader to 
follow through the analyses presented herein with numerical data and results. 

We see from Fig. 7 that the flight trajectory point for Mach 3.96 is well within the facility 
capability zone and therefore, in principle, this flight trajectory can be simulated exactly in the 
HFFAF ballistic range. In Fig. 8, we show the flight trajectory again, along with the facility 
capability zones for Mach 2.37 and 3.96 for a 3.56 cm diameter maximum model size. We see 
that the flight trajectory point at Mach 2.37 is again well within the facility capability zone. In 
principle, all of the flight trajectory points can be simulated exactly in the HFFAF ballistic range. 

This, however, is not feasible in practice because it would necessitate different model mean 
densities, diameters and sabot designs for each Mach number, resulting in an extremely 
expensive and time-consuming test entry. In the present example, since the critical event, 
parachute deployment, occurs at the lower end of the Mach number range, it was decided to try 
to provide a close range match to the flight condition at Mach 2.37. For this condition, the 
calculated 2/D for pitch oscillations is 331. Taking 4 to be 11.43 m (giving two oscillation cycles 
in the range), D would be 3.45 cm. We chose a slightly larger value of D, 3.56 cm, which is 
80% of the gun bore diameter. (The larger the model diameter, the less likely the model is to be 
disturbed or kicked by the enclosing sabot when the sabot is discarded.) With this model 
diameter, matching the Reynolds number gives a required pressure in the range of 170 Torr. 
Matching the density ratio gives a model density of 7.35 g/cm*. Hence, steel, with a density of 
7.80 g/cm’, was chosen as the model material. This defines one point on the green curve in Fig. 
8. The remaining points on the green curve were defined by choosing the range pressure for 
each Mach number to keep the same pitch oscillation wavelength — that is, to maintain 2 pitch 
oscillations in the length of the range. By moving from the red curve to the green curve (i.e., 
from an impractical test entry to a practical test entry) one gives up the possibility of obtaining 
(in principle) an exact match of p,/p and Re between the model and the flight vehicle for each 
Mach number. 

The above procedure provides very good simulation around Mach 2.37, with parameter 
mismatches of the order of 10%. At the higher Mach numbers, the mismatches can be as large 
as 20% for p,/p and 35% for Re. We note that the actual ballistic range Mach numbers, 2.08 to 
3.65, are slightly lower than the flight trajectory Mach numbers originally studied. Figure 9 
shows an expanded version of Fig. 8, showing the variable (impracticable) model and final 
(single model) trajectories in more detail. Taking the green data point at Mach 3.0 as an 
example, it can be seen that by reducing the range pressure (moving along the purple curve), 
we can get much closer to the desired p,/p and Re values. While this may appear to be 
desirable, it reduces the number of oscillation cycles in the range and, hence, was rejected as 
an option. 

The analysis discussed above for the LDSD model with the inflatable deployed configuration 
was repeated for the LDSD model with the inflatable stowed configuration. The vehicle and the 
ballistic range model are shown in Figs. 10 and 11, respectively. The resulting trajectory point at 
Mach 3.45 and the facility capability zone for this case are shown in Fig. 12. It is seen that the 
flight data point is outside the capability zone of the facility. For this case, with A/D for pitch 
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oscillations of 661 and two pitch oscillations in the length of the range, the model diameter 
would be 1.73 cm. The required model density is 23.8 g/cm”, which is unattainable. Since a 20 
mm gun was available, a model diameter of 85% of the gun caliber, or 1.698 cm, was selected. 
This is only slightly less than the value for the flight data point in Fig. 12. The model was 
designed to be fabricated primarily of a machinable tungsten alloy with a density of 17 g/cm® 
with a steel “anti-weight” in the nose to properly locate the center of gravity. This resulted in a 
mean model density of 14.24 g/cm’. To have two cycles in the range, the required range 
pressure was calculated to be 195 Torr, which determines the range model Reynolds number 
and the range model density ratio, p,/p. These values then determine the location of the upper 
green data point in Fig. 12. The lower green data point represents the range condition for some 
shots which were made at Mach 3.15 rather than at Mach 3.45. For the data point at Mach 3.45, 
the mismatches are ~20% for density ratio and 30% for Reynolds number. 


6. Calculating range simulation capabilities: example 2 — the High- 
Energy Atmospheric Reentry Test (HEART) 


The analysis of Sec. 5 was repeated for the proposed HEART vehicle. The HEART vehicle is 
shown in Fig. 13 and is a low density decelerator comprised of a centerbody surrounded by an 
inflatable cone made up of 12 tori. As will be seen, with the HEART vehicle and trajectory, there 
is more difficulty in obtaining dynamic similarity and Reynolds number matching than with the 
LDSD vehicle and trajectory. The normalized masses (mass divided by diameter cubed) of the 
HEART vehicle and the LDSD vehicle with the inflatable deployed are both in the range of 6 to 8 
kg/m°. Thus, the vehicles are rather similar to each other and the difference in the difficulty in 
obtaining dynamic similarity and Reynolds number matching is primarily due to the different 
trajectories for the two vehicles. 

The HEART flight test will have as primary objectives testing the robustness of the fabric tori 
with respect to heating and aerodynamic forces. The static and dynamic stability of the vehicle 
are also of importance and can be assessed in the ballistic range prior to the flight test. The 
Mach numbers of interest for the HEART ballistic range stability tests are 0.8 to 2.0. Figure 14 
shows the flight trajectory of the HEART vehicle and the simulation capabilities of the HFFAF 
range. Figure 14 is analogous to Fig. 7 for the LDSD vehicle with the inflatable deployed, except 
that facility capability data for the two extreme Mach numbers are shown in Fig. 14. We see 
that, unlike the case for the LDSD vehicle with the inflatable deployed, the trajectory of the 
HEART vehicle lies well outside the capability of the HFFAF range, especially if one is limited to 
using the 44 mm gun (the Ames powder gun with the most extensive user database). The 
HFFAF range, particularly with the 44 mm gun, falls well short of the flight Reynolds number 
requirements and far off the flight p./o values. With a decelerator flying at a low p,/p, the 
oscillations in pitch are relatively much more rapid than with a decelerator flying at a high p./p, 
and the A/D values are correspondingly much smaller, leading to the problem that there are far 
too few data stations per cycle to capture the pitch oscillations accurately. A measure of the 
capability of a range to capture the pitch oscillations is the station spacing divided by the largest 
model diameter (which we have been taking as 0.80 times the caliber of the largest gun). 
Characteristics of the Ames HFFAF range and three 1 atm open ranges are given in Table 1 
below. (The Ames GDF range is located in the same complex as the Ames HFFAF range; the 
Eglin ARF range is located at Eglin Air Force Base, Florida and the DRDC range is located at 
DRDC, outside of Quebec City, Quebec, Canada. As of November, 2016, the Eglin range is no 
longer operational.) 
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Table 1. Characteristics of various ballistic ranges. 


Max model Ax/D 
diameter, D 
(cm) 


Range Range length | Station spacing, 
Ames HFFAF 
Ames GDF 0.098" 


Ames GDF 0.098" 
Eglin ARF 3.38 


(m) (m) 
AmesHFFAF | 22.9 | 1.52, | 44 | 8H | 4 
[eet 
a 
| 26, | 


un 

m) 
44 
61 
44 
61 
76 


DRDC 
DRDC 12.40 


*Small station spacing achieved by having each of 8 cameras take 12 — 14 images 
**155 mm gun is currently rifled, would have to be smooth bored to do aerodynamic range work 
tAverage station spacing for first 91 m of range length 


G 
m 
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In regard to capturing short wavelength pitch oscillations, the ARF range has less capability 
than the Ames HFFAF, while the DRDC range has some advantage over the HFFAF range, 
particularly if the 155 mm gun could be brought into service. In this respect, the Ames GDF 
range has a very great advantage over the Ames HFFAF range. The Ames GDF range does 
have limitations due to its short length, however. An important comparison will be to place the 
capabilities of the other three ranges on the graph of Fig. 14. However, before doing that, we 
discuss how the capabilities of a 1 atm range appear on the p,/p — Re plane. 

Figure 15 shows the capabilities of the Ames GDF range for the HEART vehicle at the Mach 
0.8 flight data point, using the Ames 44 mm powder gun. These capabilities were calculated for 
a single material ballistic range model, assumed for model robustness. The limiting lines are 
established, in part, by limits on model diameters and pitch oscillation wavelengths. The 
maximum model diameter is taken to be 3.56 cm, which is 80% of the gun diameter. The 
minimum model diameter is taken to be 1.27 cm, based on having sufficient resolution in the 
photos of the model. (For the GDF ICCD camera system, this corresponds to having 10 pixels 
across the diameter of the model, which was judged to be the minimum to obtain acceptable 
position and angle data quality.) The maximum wavelength is taken to allow 2 cycles within 
length of the range. The minimum wavelength is taken to be that which allows data points at 6 
stations per cycle. 

The horizontal boundaries in Fig. 15 are established by the maximum and minimum model 
diameters, with the known model velocity, air density and viscosity. One criterion for the left side 
vertical boundary is established using the minimum wavelength and the maximum model 
diameter and then translating to p;/p using Eq. (3). Another criterion is p./p based on the 
assumed lowest model density of 0.94 g/cm® (polyethylene). The larger of these two values is 
used to establish the left side vertical boundary. One criterion for the right side vertical boundary 
is established using the maximum wavelength and the minimum model diameter and then 
translating to p,/p using Eq. (3). Another criterion is p,/p based on the assumed highest model 
density of 13.0 g/cm® (machinable tungsten alloy with sufficient margin for an anti-weight to 
position the center of gravity in the correct location). Whichever of these two values is the 
smaller is used to establish the right side vertical boundary. 

For the inclined boundaries, we take Eq. (7), cube both sides of the equation and then divide 
both sides of the equation by (p,/p)*. This yields 


p, _{ 1 AMpc (10) 
p \K Reu 
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Note that p is a constant for a 1 atm range. Also, note that p./p now varies as (1/Re)*, whereas 
in Eq. (7), ps/p varies as (1/Re)*”. For the upper inclined boundary, Eq. (10) is evaluated at the 
maximum A and at the lower inclined boundary, Eq. (10) is evaluated at the minimum i. 

Figure 16 shows the HEART flight trajectory (red) and the simulation capability of the Ames 
HFFAF range for the Mach 2 flight condition in blue. The simulation capabilities of the Ames 
GDF range are shown in green and those of the DRDC range in purple. Note that the simulation 
capability box of the DRDC range using its 155 mm gun does not show 6 boundary lines like the 
capability box shown in Fig. 14, since the relative positions of the limit lines causes some of the 
them to disappear. Also, for the DRDC range using its 110 mm gun, the lines squeeze together 
to such an extent that essentially only one capability point remains. Because of the high Ax/D 
value for the Eglin ARF range (see Table 1), no simulation capability of that range appears in 
Fig. 16. A relaxation of one of our criteria (such as allowing less than 6 data points per pitch 
oscillation cycle) would be necessary for a simulation capability of the ARF range to appear in 
Fig. 16. Figure 17 shows the corresponding plot for the Mach 0.8 flight condition. 

For the Mach 2 flight condition, we consider the best facility operating condition within the 
simulation capability of each range/gun combination. From Fig. 16, we then note a progressive 
improvement in the Reynolds number and density ratio simulation as one moves up in gun size 
from 44 mm to 61 mm to 110 mm to 155 mm. Within the HFFAF capabilities, the 61 mm gun 
does substantially better than the 44 mm gun. The DRDC range can come quite close to the 
flight conditions. The Ames GDF range can match the density ratio but would be a factor of ~4 
short of the flight Reynolds number. More details relating to the preceding discussion are given 
in Table 2. In Tables 2 and 3 below, R = p,/p, ps, is the ballistic range model density and P, is 
the range pressure. (R is used for clarity in Tables 2 and 3.) 


Table 2. Key parameters for simulation of HEART flight condition at Mach 2.0 in various 
facilities with various guns. 


(mm) 
HFFAF 61 


3.49x10° 7270 1.74x10° 13,800 
DRDC 3.49x10° 7270 4.07x10° 11,340 
DRDC 155 3.49x10° 7270 5.00x10° 7270 


44/20 3.49x10° 7270 0.83x10° 7270 


Table 2 (continued). Key parameters for simulation of HEART flight condition at Mach 2.0 


in various facilities with various guns. 
Range Gun D, Osr P, Mismatch Mismatch 
Tele lool eS 
(mm) m) (g/cm?) (Torr) Re,/Re, R/Ri 
HFFAF 61 
8 


( 


3.56 13.0 319 9.18 3.57 
4.88 13.0 600 2.00 1.90 


C 
DRDC 110 8 13.5" 760 0.858 1.56 
DRDC 155 11.2 8.68 760 0.698 1.00 
44/20 1.85t 


tCould perhaps use 20 mm gun 
*If ps = 13 g/cm’, only 5.88 points per oscillation cycle would be obtained 
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For the Mach 0.8 flight condition, as for the Mach 2 flight condition, a progressive 
improvement in Reynolds number and density ratio simulation as one moves through the 
increasing diameters of the Ames guns in the HFFAF range to the DRDC guns in the DRDC 
range is noted. However, in this case, even the best DRDC simulation is short of the flight 
Reynolds number by a factor of 2.2 and off density ratio by a factor of 3.3. The Ames GDF can 
match the flight density ratio, but is well short of the flight Reynolds number (by a factor of 7.6). 
More details relating to the preceding discussion are given in Table 3. 


Table 3. Key parameters for simulation of HEART flight condition at Mach 0.8 in various 
facilities with various guns. 


Range G Rey Ry Re, R,; 
(mm) 


DRDC 110 5.04x10° 2034 1.63x10° 13,230 
DRDC 155 5.04x10° 2034 2.29x10° 6690 
| GDF | 44¢ ~~ | 5.04x10° 2034 0.66x10° 2034 


Table 3 (continued). Key parameters for simulation of HEART flight condition at Mach 0.8 


in various facilities with various guns. 
Range Gun D, Osr P, Mismatch Mismatch 
Te lalele =e 
(mm) (cm) (g/cm?*) (Torr) Re;/Re, RJR» 
HFFAF 61 
8.8 


un 
HFFAF 44 5.04x10° 2034 0.231x10° 30,260 
HFFAF 61 5.04x10° 2034 0.596x10° 16,090 
44t 


4.88 13.0 014 8.45 7.91 
DRDC 155 12.4 


15.8" 760 3.10 6.50 
7.99 760 2.20 3.29 


*If 0; = 13 g/cm’, only 5.44 points per oscillation cycle would be obtained; if ps = 14.24 g/cm° 
could be obtained, as was done for the LDSD stowed model, then the number of points per 
cycle would increase to 5.70. 


Tables 2 and 3 give the required size of the model (D,) for each possible condition of choice. 
Also given are the required model densities. Models with densities of 13 g/cm® (and above) 
would likely be fabricated mainly out of machinable tungsten alloy, with a density of 17 g/cm’. 
Several ballistic range models in the tables have required densities of 8 — 8.7 g/cm’. These 
could be fabricated mainly of steel (7.8 g/cm’*) or brass (8.7 g/cm’). In one case, the required 
model density is 2.43; this model could be fabricated mainly of aluminum (2.7 g/cm’). Internal 
hollows in the model or magnesium (1.8 g/cm”) inserts could bring the density down to that 
required value. 

The choice of which range (and gun) to be used will depend upon many factors, two of which 
are: 


Relative importance placed on Reynolds number and density ratio mismatch. (The 
density ratio mismatch translates to a mismatch in A/D, the normalized pitch 
oscillation wavelength.) 


The increasing expenses for the models, sabots and gun operation as the gun size 
increases. 
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The best choice of range and gun for HEART type vehicles depends upon the relative weighting 
given to the various factors and hence, we make no attempt to give a definitive answer here. 


7. Determination of decrease in velocity and Mach number per cycle 
of pitch/yaw oscillations and during traverse of the range test section 


If the model has the same r,/D as the flight vehicle, and 4/D is matched, the fractional velocity 
drop per pitch oscillation cycle will be the same for the flight vehicle and the model. Thus, the 
deceleration of the flight vehicle and the model will be dynamically similar. For some vehicles, 
particularly those with inflatable or deployable aerobrakes, it may not be feasible, with a robust 
gun-launched model, to have r,/D the same for the flight vehicle and the model. The fractional 
velocity drop per cycle can be calculated from the following equations, starting with Newton’s 
law, force (F) = mass (m) x acceleration (a). The aerodynamic drag force is given by 


] 
F=CpA- pu’ (11) 
The mass times acceleration term is given by 
du du 
ma = m— = mu— (12) 
at dx 


where t is time and x is distance along the flight trajectory. From Eqs. (11) and (12), we get 


du _C,Ap 


dx = ae 


dx (13) 
u 2m 2V p, 


Evaluating du/u for dx = the wavelength, A, using Eq. (2), we get 


d AD 2 91 
a ee {Oe ete eee (14) 
u |, V p,C,, D 
Applying Eq. (2) once more, we get 
2 
du he 5. a D (15) 
u |, C, \D]} A 


For the LDSD vehicle with the inflatable deployed and the conditions discussed in Sec. 5 with 
Mach numbers from 2.37 to 3.96, the fractional velocity drop per cycle is ~0.04 for the flight 
vehicle and ~0.05 for the ballistic range model. These differences between these small velocity 
drops very likely has a negligible effect on the simulation of dynamic similarity between the 
model and the flight vehicle. For the HEART vehicle, with a large inflatable, the velocity drops 
per cycle are larger, being ~0.14 for the flight vehicle and ~0.20 for the model. These 
differences still may be acceptable for sufficiently accurate dynamic similarity. 

The fractional velocity drop may be reduced by decreasing the model diameter D, while 
maintaining 2. However, this moves one away from A/D pitch oscillation similarity and reduces 
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the Reynolds number, both of which are undesirable. For example, let D be reduced to 0.8 of its 
original value, while 4 is maintained, which can be done by decreasing p_ to 0.64 of its original 
value. This would reduce the fractional velocity drop by 20%, but would also put A/D off dynamic 
similarity by 20%. Further, the Reynolds number would be reduced to 0.51 of its original value, a 
49% reduction. 

If it is desired to move the model r,/D closer to the flight vehicle r,/D to achieve closer 
dynamic similarity of the velocity drop, this can, in principle, be done by building a two-material 
model, say of magnesium and steel. Such a model would be much more complex than a single- 
material model and may be difficult to make sufficiently robust to survive the gun launch 
process. There may be a strong tendency for such a model to come apart at the joint between 
the two materials. 

Reference 3 gives other useful velocity drop criteria. One criterion of this reference is that the 
Mach number change (drop) on traversing the range test section should be less than 0.2. We 
take the fractional velocity drop per oscillation cycle from Egn. (15) and transform to get 


dM 


Fr = (16) 


range 
where M is the Mach number and L is the length of the range. From Eq. (16), we get 


du 


DL 
u 


dM ad 
AD 


(17) 


range range 


For the LDSD model with the inflatable deployed, for the conditions discussed in Sec. 5 with 
Mach numbers from 2.37 to 3.96, the Mach number drop on traversing the range from Eq. (17) 
ranges from 0.24 to 0.32. These drops somewhat exceed the suggested limit of 0.20 of Ref. 3. 
For the HEART model shown in Table 2 of Sec. 6, in the NASA Ames GDF facility over the 
Mach number range of interest (Mach 2.0 to Mach 0.8), the Mach number drops on traversing 
the range are much larger, ranging from 0.52 to 0.58. This model is a relatively light model and 
slows down very rapidly — hence, the large Mach number drops are unavoidable, unless one 
moves well away from dynamic (A/D) similarity between the model and the flight vehicle. 

A second criterion given in Ref. 3 is that the Ax due to deceleration in the range, i. e., the 
difference between the distance travelled by the actual model and the distance travelled by a 
model with zero deceleration, be greater than a certain value, which is given as 2.5 cm for the 
NASA Ames HFFAF facility. This limit is a lower limit for measurability. This limit corresponds to 
accuracy of 0.013 cm in measurements of the model position in the shadowgraphs and a time 
measurement accuracy of 0.02 microseconds. Ax can be calculated from Eq. (18): 


L (18) 


range 


For the LDSD model with the inflatable deployed and the HEART model discussed in the 
preceding paragraph, Ax is between 60 and 160 cm, much larger than the minimum Ax criterion 
(2.5 cm) given in Ref. 3. Hence, for these models, there will be no problem regarding this 
criterion. 

For the LDSD vehicle with the inflatable deployed and the range simulation of same, the key 


du 
parameters necessary for the calculation of — 
u 


, AMrange and Ax and the results themselves are 
A 


given in Appendix B. 
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8. Determine swerve due to oscillating lift caused by pitch (or yaw) 
oscillations 


Following the treatment of Ref. 4, we estimate the swerve due to the oscillating lift caused by 


pitch (or yaw) oscillations. This swerve will be compared to the estimated precision with which 
the model position can be measured to assess the accuracy with which the lift coefficient slope ( 


Ci ) can be determined. We consider the pitch angle of attack (a) to be oscillating sinusoidally: 


a=. sin( 7) (19) 
A 


where x is the distance along the range. The lift (L) at an angle of attack a is given by: 
Ll 
Ls ae A A (20) 


Combining Eqs. (19) and (20), we get: 

_(2nx\1 > 
L=C, a... sin} ——|—pu'A 21 
La ~ “max | X Se ( ) 


Using Force = Mass x Acceleration, we get 
y=L/m (22) 


where y is the acceleration perpendicular to the flight path. Combining Eqs. (21) and (22) and 
using x = ut (assuming constant velocity), where t is time, we get 


y . pu’ A (23) 


oe Cy, Bmax . 2et) 1 
= — — Sl 
Mm 


Integrating twice and evaluating at the maximum of the sine function, we get 


-_ C, Ob max A*pA 


D nite = 2 


(24) 


m S70 


where Ymax IS the maximum swerve due to the oscillating lift. Substituting for A from Eq. (1), we 
get: 


a (25) 


Using J = mr. , we get 
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Cc 2 
a a |B] (26) 


Me 


Using Eq. (26), we can evaluate Ymax for the LDSD model with the inflatable deployed and the 
HEART model discussed in the previous section. AS was mentioned previously, the accuracy of 
the position measurements in the NASA Ames HFFAF facility is given in Ref. 4 as 0.013 cm. 
For an example given in Ref. 4, with an max of 10 degrees, a Ymax Of 0.5 cm is calculated. With 
the measurement accuracy of 0.013 cm, the maximum swerve displacement was about 40 
times the measurement accuracy. For the LDSD model in the NASA Ames HFFAF facility, Ymax 
was estimated as 0.19 to 0.25 cm, being 15 to 19 times the measurement accuracy. For the 
HEART model in the NASA Ames GDF facility, ymax was estimated to be 0.13 to 0.18 cm. If the 
measurement accuracy of the Ames GDF facility is taken to be equal to that for the Ames 
HFFAF facility (0.013 cm), the estimated ymax values are 10 to 14 times the measurement 
accuracy. These ymax values are not as favorable as those in the example given in Ref. 4 and 
point to the difficulty in getting good lift curve data. The problem becomes more severe at 
smaller amax values, such as 5 degrees or 2 degrees. Unfortunately, the parameters Cs 


Cn, and r,/D in Eq. (26) are characteristic of the model and cannot really be changed 


significantly. D is usually constrained by the size of the launchers and by the need to have an 
adequate number of oscillation cycles in the range. 

For the LDSD vehicle with the inflatable deployed and the range simulation of same, the key 
parameters necessary for the calculation of Ymax and the ymax values themselves are given in 
Appendix C. 

It is concluded that for testing the LDSD model with the inflatable deployed in the HFFAF and 
the HEART model in the GDF, errors of the order of 5 to 10% in the lift coefficient slope may 
occur due to measurements inaccuracies when Gmax Is ~10 degrees. Errors would be greater at 
smaller Gmax Values and less at larger Omax values. In such cases, one needs to make every 
effort to get accurate swerve measurements. 


9. How to determine swerve due to lift for lifting models 


If the model lifts, one can determine the (vertical) swerve due to lift down the length of the 
range. The equations to determine the swerve follow. Again starting with Newton’s law F = ma, 
the aerodynamic lift force is 


1 
F= CAS pu’ (27) 
The lateral acceleration a is, then 
1, 
C,A— pu 
q=—__4_ (28) 
m 


The swerve, y, over a time of flight t over a range length L is, assuming constant u, 


i 
2 ou 


Il. 
=—at = 29 
y 5 (29) 
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Combining Eqs. (28) and (29) and using m = p,V, we get 


1 LAD 
yore, =F (30) 
4 “DV p, 
Using Eq. (2) with Eq. (30), we obtain 
2 2 
y oC, dD I 
EC DAD oo 
and 
r DL 
702 OED = 


Note that the length, L, in equations (29) to (32), must be the length from the gun muzzle to the 
last data station rather than the distance from the first to the last data station. For the Ames 
HFFAF, L is 32.92 m. For the Ames HFFAF, the maximum lift swerve allowable before the 
model lifts out of the field of view of the last shadowgraph station is 17 cm. By tilting the gun 
downwards (at angles up to ~1 degree) so that the model first dives down, then levels out and 
then starts to rise upwards along the length of the range, the maximum allowable lift swerve can 
be increased by almost a factor of three to 47 cm. In three test campaigns®®, by tilting the gun, 
models with lift coefficients of 0.3 - 0.6, ~0.35 and ~0.4 were successfully tested in the Ames 
HFFAF. The corresponding lift/drag ratios ranged from 0.25 to 0.51. It should be pointed out 
that the tests of Refs. 6 and 8 were performed in the Ames HFFAF, while the tests of Ref. 7 
were performed in the Ames GDF facility (see Sec 6). At the present time, it is only possible to 
tilt the powder guns and not the two-stage light gas guns. 

Models with higher L/Ds may be tested by decreasing the model diameter D, decreasing p or 
increasing ps; (or a combination of these steps), but this moves one away from A/D pitch 
oscillation similarity and reduces the Reynolds number. Due to the small size of the ballistic 
range models, the Reynolds numbers are often already somewhat below the flight values. 
Hence, a further reduction in Reynolds number is usually undesirable. Table 4 shows the 
normalized variables for several techniques which could be used to reduce swerve (s = 
normalized swerve). Cases 1 to 4 show techniques for reducing the swerve by 20%. For case 1, 


Table 4. Normalized variables (baseline variables = 1) for several techniques to reduce 
swerve. 


Rael a Pl 


|Caset | 1 | 1 | 4 ~+| 08 | 08 | 14 | 125 | 08 | 08 


064 | 1 | 141 | 125 | 1.25 | 08 | 064 | 08 
|Case3 | 156 | 1 | 156 | 141 | 1.25 | 1.25 | 08 | 1 | 08 


the diameter is reduced by 20%, which leaves the dynamic similarity unchanged, reduces Re by 
20% and increases the number of cycles by 25%, which may or may not be acceptable. For 
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case 2, the same reduction in the swerve is achieved by decreasing the gas density. Here the 
penalties are more severe, A/D is increased by 25%, Re is decreased by 36% and the number 
of cycles is reduced by 20%. The reduction in the number of cycles may result in less than 4 
extrema of the pitch oscillation angles within the length of the range and a consequent reduction 
in the accuracy of the determination of the growth or decay rate of the pitch angle oscillations. If 
it is possible to increase ps by 56%, one may arrive at Case 3. Here, Re is unchanged, but A/D 
increases by 25%, which decreases the number of cycles by 20%. For case 4, both p and D are 
decreased, which permits one to maintain 4 and the number of cycles. The penalties are a 28% 
reduction in Re and a 12% increase in A/D. Since the cycles can be maintained for the 
technique of case 4, that case was extended to achieve a more substantial reduction in s (by a 
factor of 2) in case 5. Here, the Reynolds number penalty is severe, by almost a factor of 3. 
Also, A/D is increased above its baseline value by 41%. 

Models with much higher L/D values can be tested in the Ames GDF facility, but this range is 
a short (4.88 m) open (1 atm) ballistic range, which limits model size and the ability to achieve 
dynamic and Reynolds number similarity. 

One method of making a model lift is to insert a weight (and/or “antiweight” — a weight with a 
density less than that of the bulk model material) which is located off the centerline of the model. 
This causes the model to trim at an angle of attack, producing lift. In some cases, when the 
weight was placed above the centerline of the model, the model was found to roll and, hence, to 
lift sideways out of the field of view. This problem was solved by reversing the model so that the 
weight was below the model centerline. This change, of course, reversed the direction of lift of 
the model. 

For the LDSD vehicle with the inflatable deployed and the range simulation of same, the key 
parameters necessary for the calculation of y and the y values themselves are given in 
Appendix D. To illustrate this procedure, in the Appendix, we have assumed that C, = 0.23. The 
actual C, for this vehicle is zero and hence y for the actual, non-lifting vehicle (and model) would 
be zero. The vehicle (and model) would have to be weighted off-center to produce a non-zero 
CL. 


10. Model design — materials, weights 


Table 5 shows materials frequently used to fabricate ballistic range models. Most of these 
materials have been successfully used by the present author. In most cases, a material can be 


Table 5. Materials for ballistic range models 


Material Compound, alloy, remarks 
g/cm? 


(UHMW = ultra high molecular weight) 
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chosen which has a density close to the desired density calculated in sections 5 and 6. In some 
cases, two-material models, screwed or press fit together may be required to achieve the 
desired average model density. 

In many cases, a two- (or three-) material model is required to properly place the center of 
gravity (CG) of the model in the correct axial location. To do so, a nose weight or “anti-weight” is 
often required. The nose weight may be press-fit into the main model body or the model can be 
made in several pieces and screwed together. Figure 18 (see Ref. 5) shows a model with a 
press-fit nose “anti-weight”. Similarly, to offset the CG from the axial centerline so as to 
generate lift, weight/“anti-weights” can be press-fit sideways into the main body of the model. 
Figure 19 (see Ref. 6) shows a model with such a press-fit weight and an “anti-weight” to 
laterally offset the CG. The same type of CG offset can be obtained by a sideways offset weight 
trapped by a threaded retainer as shown in Fig. 20 (see Ref. 7). We note that by proper 
selection of the axial positions of the weights and “anti-weights” shown in Figs. 19 and 20, the 
presence of these weight/“anti-weights” can also be used to adjust the axial location of the CG 
as required. Figure 21 (see Ref. 8) shows a three piece tabbed model. Several methods of 
adjusting the CG location in this model are possible. Small changes in the axial position of the 
CG can be made by adjusting the sizes of the cavities A and B. A large change in the axial 
position of the CG was made by changing from a brass forebody to an aluminum forebody while 
also changing from a ABS plastic afterbody to an aluminum afterbody. The CG is obviously 
offset laterally by the presence of the tab. However, considerable additional lateral CG offsets 
can be made by adjusting the sizes of the cutouts C and D. 


11. Sabot design — materials 


Table 6 shows materials used for sabots for ballistic range models. The first three materials 
have frequently been used and have been successfully used by the present author. The ST-801 


Table 6. Materials for ballistic range sabots 


Material Compound, alloy, remarks Density, 
g/cm? 


ABS, high impact Acrylonitrile-Butadiene-Styrene 


ST-801 Nylon Poly-amide, DuPont 
Polycarbonate 
Radel R5500 Polyphenylsulfone 


“super tough” Nylon has, in general, given us good service. However, unexpected failure has 
occasionally occurred with this material. The impact strength of plastics are frequently 
measured by the notched Izod impact test (ASTM D256). The Izod value for ST-801 Nylon is 
quoted by DuPont as 17 — 20 ft-Ibf/in. Izod tests that were made for doubtful batches of Nylon 
have given values as low as 5 — 6 ft-lbf/in. Thus, for severe launch conditions, it may be 
desirable to have Izod tests performed before fabricating components from ST-801 Nylon so as 
to ensure their performance. As a standard practice in the HFFAF, at least 5 Izod specimens 
are made and tested. Ideally, some specimens should be cut in all three possible orientations. 
However, it is realized that for round stock or flat stock with a maximum dimension of less than 
2.5°, only two orientations of the Izod specimens are possible. 

For moderate launch conditions, ABS has given us good service. One must be sure to obtain 
“high impact” ABS with Izod values of 5 — 7 ft-lbf/in, rather than “medium impact” ABS with Izod 
values of 2 — 4 ft-lbf/in. The tensile strength of ABS is given as about 6 ksi, nearly the same as 
that of ST-801 at 2.3% water content. (This water content occurs when the Nylon is in 
equilibrium at 50% relative humidity air.) The density of ABS is slightly less than that of ST-801. 
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Our current policy is to use ABS (and sometimes ST-801) for moderate launches and ST-801 
for more severe launches. If the ST-801 is tested and proven to have the claimed Izod values, it 
has roughly 3 times the impact strength of ABS and thus, is the sabot material of choice for 
severe launch conditions. The Radel has successfully been used at the NASA Marshall center, 
but not by the present author. It has a tensile yield strength of 10 ksi and an Izod _ strength of 13 
ft-lbf/in, but a somewhat higher density than the other plastics in Table 6. 

Many designs for sabots made of polycarbonate (Lexan, Zelux) are given in Ref. 9. 
Reference 10 states that “Polycarbonate plastic (Lexan) is widely used because it has a higher 
impact strength and a high(er) modulus of elasticity than most other plastics.” Our experience 
with polycarbonate is that it can easily craze crack if it is not carefully double annealed. 
Polycarbonate is also subject to degradation due to age, ultraviolet light and some lathe cutting 
fluids. We are not currently using polycarbonate for sabots, although its nominal strength is 
greater than that of ABS or ST-801 Nylon. If one chooses to use polycarbonate for sabots (or 
models) it would be prudent to use new material that has been double annealed. 

The mass of the sabot is sometimes called “parasite mass” because the goal is to launch the 
model, but, in most cases, a sabot is necessary. The sabot can make a substantial contribution 
to the total launch package mass. To enable higher velocities to be obtained or to reduce the 
stress on the gun and the launch package, a lighter sabot would be preferred. Because ABS 
and ST-801 Nylon have the lowest densities in Table 6, the sabot parasite mass will be less 
when the two former plastics are used, instead of polycarbonate or polyphenylsulfone. 


12. Sabot separation 


A typical launch package consists of a model encased within four (or two) interlocking sabot 
“fingers”. Figure 22 shows a model inside a four-fingered sabot with one of the fingers off to one 
side. When the launch package exits the gun muzzle, the sabot fingers must be separated 
(move apart sideways) from the model, allowing it to fly freely. Ideally, the model then passes 
through a center aperture in a steel cone (the “sabot stripper”), while the sabot fingers impact 
the outer surface of the cone and are stopped. Figure 23 shows the arrangement of the gun, 
dump tank, sabot stripper and shadowgraph test section for the NASA Ames HFFAF facility. 
(The use of the “slap” paper will be discussed in Sec. 14.) The sabot fingers can be driven apart 
sideways by several techniques. One technique consists of shaping the front of the sabot 
fingers to form a conical (typically at a half angle of ~45 degrees) ramp-like surfaces. The 
aerodynamic ram pressure then will open up the sabot fingers from the front. This is called 
aerodynamic sabot separation and can be seen in the sabot design shown in Fig. 22. As a 
rough rule for the Ames HFFAF ballistic range, a minimum pressure of 30 — 50 Torr of air is 
required for adequate aerodynamic sabot separation. If a cylindrical hole is cut in the rear part of 
the sabot fingers, the muzzle blast of the gun can be used to open the sabot fingers from the 
rear. This is called muzzle blast sabot separation. At times, both separation techniques can be 
used together. Figure 24 shows a launch package using both types of sabot separation 
techniques. Figure 25 shows several progressive frames from a sabot separation video for the 
launch package shown in Fig. 24. Figure 26 shows frames from a sabot separation video for a 
launch package identical to that shown in Fig. 24 except without the muzzle blast cavity. It is 
evident that with the muzzle blast cavity, the sabot finger separation angles (from the rear) are 
much larger than those without the cavity. 

When launching into very low (e.g., < 30 - 50 Torr) test section pressures, as is typically done 
for calibration spheres, the aerodynamic separation technique is not effective. Hence, muzzle 
blast separation is the technique of choice. (In some cases, while the initial launch is into a 
vacuum, farther down range the model passes through a diaphragm into the pressurized region 
of the range, where the aerodynamic studies take place.) Figure 27 shows a launch package 
designed for launching a very light model into a vacuum using muzzle blast sabot separation. 
Although the forward ends of the sabot fingers are cut into a 45 degree cone shape, there are 
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no aerodynamic separation forces since the launch is into vacuum. In addition, aerodynamic 
separation could not be used in this case, because the aerodynamic forces would tend to drive 
the very light model back into contact with the sabot fingers, resulting in unacceptable wild 
motions of the model. 

In some cases, it is advantageous to place an obturator (gas seal) cup behind the sabot 
fingers to block the passage of drive gas up through the joints between the serrated sabot 
fingers. This prevents the drive gas from heating the model for aerothermal tests. However, this 
prevents the use of muzzle blast sabot separation. The obturator cup generally follows the 
model down the range, but drops behind sufficiently so as not to affect the flight of the model. 
Figure 28 shows a launch package with an obturator cup (beneath the sabot fingers) designed 
for aerodynamic separation of the sabot fingers. Discussion of obturator cup design and tests is 
presented in Ref. 11. 

In some ranges, (e.g., the Ames Vertical Gun Range - AVGR), the sabot fingers are 
separated by centrifugal forces, using a rifled gun barrel. This will not work for the Ames HFFAF 
and GDF ranges, since the guns used in these ranges are not rifled and since rapidly spinning 
models are not suitable for vehicle aerodynamic studies, wherein the flight vehicles are not 
spinning. 

Reference 12 has a good discussion of sabot separation techniques and mechanical sabot 
separation techniques in particular. One mechanical sabot separation technique is to have a 
reduced (less than the barrel diameter) diameter hole plate at the muzzle to intercept the outer 
edge of the sabot and slow it down or stop it. In one version, the hole plate successively 
stepped down to 82% of the barrel diameter and stopped the plastic sabot. This was used up to 
launch velocities of ~200 m/sec. (Note that this technique does not appear in Ref. 12.) A 
second version, which does appear in Ref. 12 (with details in Ref. 13) has a hole plate of 90% of 
the barrel diameter which intercepts the outer part of a foam sabot enough to slow it down to 
drop behind the model. This was used with launches at velocities up to 800 m/sec. Reference 
12 also discusses contracting vented barrels and deflecting barrels and ramps. 


13. Stress calculations 
13.1. Launch stresses 


As a first step, one needs to calculate the muzzle velocity of the gun. From the maximum 
desired mid-range model velocity, using the CFD-estimated drag coefficient, one can back 
calculate to the maximum muzzle velocity required for a given model in a given test campaign. 
Then, given the length of the gun barrel, the average acceleration needed to achieve this 
muzzle velocity is calculated. The maximum acceleration in the gun is typically about 2.5 times 
the average acceleration for powder gun launches, but can be as much as 3 to 6 times the 
average acceleration for two-stage light gas gun launches. In some cases, for powder gun 
launches at NASA Ames, powder chamber pressure histories for a substantial number of 
different shot conditions are available. These data, with some scaling, can be used to estimate 
the maximum model base pressure and hence, the maximum acceleration for the shot in 
question. For two-stage light gas gun launches, it may be possible to run a CFD simulation of 
the shot in question to get a value for the maximum launch package acceleration in the barrel. 
With the maximum launch package acceleration in hand, one can proceed to the calculation of 
the model and sabot stresses. 

We first consider a conventional (compact) model in a sabot. Figure 29 shows the NASA 
Ames Mars Science Laboratory (MSL) model in the sabot®. The arrows in Fig. 30 show the high 
bearing stresses under the model, the lower (near hydrostatic) pressures in the sabot away from 
the model base and the drive gas pressures on the base of the sabot. These combined stresses 
will lead to shear stresses (see Fig. 31) tending to allow the model to “fall” through the sabot 
due to the very large launch accelerations. The difference between the bearing stresses 
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beneath the model and the near hydrostatic pressures outside the model base should be 
compared with the compressive strengths of the model and sabot materials. The shear stresses 
shown in Fig. 31 should be compared with the shear strength of the sabot material, which can 
be estimated as about 60% of the tensile strength if actual shear strength data is not available. 
The stress analyses can be performed in the older manner as described above, or a finite 
element model can be constructed (e.g., in SolidWorks) and the von Mises stresses found. The 
maximum von Mises stresses can then be compared with the uniaxial yield stresses of the 
relevant model and sabot materials. 

Figure 32 shows a ballistic range model of a vehicle with a compact centerbody and a 
relatively thin aerobrake (a solid, fixed structure designed to simulate an inflatable or 
deployable) in its launch sabot. This model simulates a flight vehicle similar to the HEART 
vehicle discussed in Sec. 6 and shown in Fig. 13. With this type of model, attention must be 
paid to additional failure modes that are unimportant with the typical compact model. 
Considering the centerbody as an isolated unit, the bearing stress under the centerbody and the 
sabot shear stress needed to prevent the centerbody from falling through the bottom of the 
sabot need to be examined. If it is concluded that the isolated centerbody cannot be supported 
by the sabot base, then the centerbody plus aerobrake considered as a single unit should be 
considered. The bearing stress for the unit as a whole can then be examined; this stress will be 
much lower and hence, more favorable than the bearing stress for the isolated centerbody. The 
sabot plastic shear stresses necessary to prevent the centerbody/aerobrake from falling through 
the sabot plastic should then be examined. Since, in this case, the aerobrake and the 
centerbody must act as a unit, the shear stress at the joint between to two components will be 
fairly high — it must be examined and determined to be within the capabilities of the material 
under consideration. If the model is comprised of, say, a magnesium aerobrake and a steel 
centerbody, the steel will tend to fall out of the magnesium and the shear stresses on the 
threads between the two materials must be examined and determined to be at an acceptable 
level for the weaker of the two materials. A finite element analysis of the model-sabot 
combination can be performed as described in the preceding discussion and the following 
stresses should be examined. 


Bearing stresses under the centerbody and the aerobrake. 

Shear stresses in the sabot material under the outer edge of the centerbody. 

Shear stresses at the joint between the centerbody and the aerobrake, particularly thread 
shear stresses if the two parts are made of different materials and threaded together. 


We note that the model-sabot combination of Fig. 32 is a concept and has not actually been 
launched to date. However, an ADEPT model with a compact centerbody and a relatively thin 
aerobrake has been successfully launched using the model-sabot combination shown in Fig. 33. 
The previous calculations are made as if the model were in a steady-state environment 
maintained at the maximum launch package acceleration. However, these extreme conditions 
are short-lived and there is considerable evidence that shows, at higher muzzle velocities (i.e., 4 
km/sec and above), at least for compact models (e.g., spheres), the launch packages can 
tolerate stresses that are nominally well above the static material yield stresses for the short 
duration of the launches. Sufficient time may not elapse during the launches for the models to 
fall through the sabots. For the launches described in Figs. 3.36, 3.39, 3.41, 3.42 and 3.43 of 
Ref. 14, and in particular for the relatively large sphere launches of Figs. 3.42 and 3.43, bearing 
and shear stresses calculated as described in the earlier paragraph regarding compact models 
are 2 to 5 times the nominal static yield stresses and yet, satisfactory launches were 
presumably obtained. The model-sabot configurations for these launch packages are shown in 
Figs. 34 — 36. Note that there is some uncertainty regarding the sphere diameter shown in Fig. 
3.43 of Ref. 14. No reference is given in Ref. 14 for this launch package and the stresses 
calculated for the larger sphere [Fig. 36 (b)] are extremely large. It might be better to key future 
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designs off the launch package with the smaller sphere [Fig. 36 (a)], where the stresses are 
more comparable to stresses for other successful designs. These launches were at velocities of 
7.3 — 9.4 km/sec. The calculated stresses for the higher velocity launches (up to 5.7 km/sec) of 
the model shown in Figs. 29 — 31 also greatly exceed the nominal static yield stresses. Such a 
favorable situation may well not be the case for lower velocity launches which last for 
considerably longer periods of time (as scaled by the launch tube diameter). Failed launches, 
where the model has fallen through the sabot rear, have been seen in X-ray photographs (e.g. 
see Fig. 37). For compact models at high velocities, it would seem desirable to key off earlier, 
successful sabot designs, even though the calculated stresses are above the static yield 
stresses of the materials. The best philosophy may be to design the sabots so that the 
calculated stresses are no greater (relative to the material yield stresses) than those calculated 
for successful sabot designs, even though both sets of stresses exceed the material static yield 
stresses. 

If stresses, calculated as described above, are found to be pushing or exceeding the limits 
(either the static yield stress limits for lower velocities or the limits of successful designs for 
higher velocities) and cannot be avoided by sabot and model design changes, then the launch 
package should be regarded as experimental. As such, a definite failure risk must be 
acknowledged, and it would be desirable to have muzzle view and separation videos and/or 
muzzle X-ray pictures to verify the quality of the launch. (It should be noted that for ANY new 
model/sabot design, it is desirable to have muzzle and sabot separation videos.) It must be 
realized that a complete failure of the sabot (i.e., model fall-through) within the launch tube 
could result in the model scraping along the gun barrel and ruining it, requiring a time- 
consuming and expensive barrel replacement. 


13.2. Aerodynamic stresses 


For models with a centerbody and a relatively thin aerobrake (see Fig. 38) in free flight, the 
aerodynamic forces will produce substantial bending and shear stresses near the joint between 
the centerbody and the aerobrake. These stresses must be calculated at the most severe 
muzzle velocity/range gas density condition for the model in question. If the model is made of 
two pieces and threaded together [(Fig. 38(b)] the joint will be able to withstand much lower 
stresses than for the case of a single piece model. If the stresses are too high, stronger 
materials or a greater aerobrake thickness near the aerobrake-centerbody joint may be 
necessary. If heavy gas jets are used to “slap” the model to induce pitch oscillations (discussed 
in Sec. 14), the aerodynamic forces must be evaluated for the heavy gas condition. 


13.3. Slap stresses 


In some cases, varying amplitude pitch oscillations are produced by “slapping” the model by 
having it impact a sheet of paper or cardstock which covers one-half of the frontal area of the 
model. (This procedure is discussed further in Sec. 14 and Ref. 5.) For bodies with thin 
aerobrakes or a thin front shell (e.g., see Fig. 39), one must also evaluate the stresses 
produced by the “slap” process. From conservation of momentum, assuming an inelastic 
collision and that the mass loading (g/cm*) of the slap media is much less than that of the 
aerobrake, the velocity of the thin aerobrake or front shell relative to the rest of the model after 
the slap impact, Ug, will be given by 


u_=—— Eq. (33) 


where 
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u = velocity of model, cm/s 

Mp = mass loading of slap paper, g/cm? 
Oa = density of thin shell, g/cm?” 

t, = thickness of thin shell, cm 


The kinetic energy per unit mass of the moving thin shell relative to the main mass of the 
model, Ka, is 


Ki = =u, (34) 


To calculate the exact shape that the thin shell will assume at maximum deflection after the 
impact is a very complex problem. We can get a rough estimate of the amount of bending 
energy in the shell by evaluating this energy for a deflection curve based on uniform pressure 
loading. We assume that this analysis gives us a curve of Mp = M,(r), where My is the bending 
moment per unit length and r is the radius out from the central axis. For the model shown in Fig. 
39(a) one might use a clamped boundary condition at the edge of the thin shell. For the model 
shown in Fig. 39(b), the boundary condition might approach simply supported to some degree 
since the thickening ring around the edge of the thin shell could tilt and pull away from the 
material of the afterbody. The shear energy per unit mass is much less than the bending energy 
per unit mass if the thickness of the shell is much less than the radial extent of the shell, varying 
as [(thickness)/(radius)]°. The maximum tensile or compressive stress at the shell surface, op, is 


My 
I, 


O,= (35) 


where y is the distance of the shell surface from the shell neutral axis, equal to t,/2 and I, is the 
moment of inertia per unit length of the shell cross section, equal to t,"/12. Hence, 


6M 
O, = 2 u (36) 


a 


The strain energy per unit volume at the shell surface, Es, is given by 


2 
O 
E =—o,€,=— 37 
7 De oe 
where €» Is the strain and E is Young’s modulus of elasticity. The average strain energy over the 
shell is one third the peak value at the shell surface. Hence, 


E,= = (38) 
Substituting for op, we get 

. 6M; 

E,= OF (39) 
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Dividing by the shell density, pa, gives the energy per unit mass 


= _ 6M, 
eae (40) 
t Eo, 
We then average over the thin shell to get 
- 6|M,| 
ee (41) 
t Ep, 
= 6 : = | 
We then evaluate E,,, at Max(o,,) = Max(M,)— = 0), , call this value E,,, ,,, and compare it 
t 
Bes yielding is possible or likely. If AK, << E vn max , one should be safe. 


with Ka. If K, = 
-E 


From K, a safety factor of perhaps 3 may apply, but this is a matter of professional 


sm,max ” 


judgment. 


For a given slap impulse density, ump, the relation of AK, with E improves linearly with 


sm,max 


decreasing E and increasing pa and improves quadratically with increasing ta and oyie. If, after 


adjusting these parameters, K,, is still too high compared with EF other options are 


available to “slap” the model. One may use a succession of thinner slap papers or jets of heavy 
or light gases to “slap” the model. For a compact model with a thin shell, one could cut a hole in 
the slap paper so that the paper impacts the model only outside the thin shell area. For compact 
models, one can avoid the thin front shell by using different model materials. Metals of choice 
are tungsten alloy, brass, steel, titanium, aluminum and magnesium. Densities of these metals 
are given in Table 5 of Section 10. Such a model design is shown in Fig. 40. For compact 
models without a thin front shell, the aerodynamic and paper slap stresses are unimportant and 
do not need to be analyzed for model design. 


14. Tilted launches and “slapping” of model to produce pitch/yaw 
oscillations 


To characterize the aerodynamics of a model fully, a number of shots must be made, with 
the model having a wide range of amplitudes of pitch oscillations. Typically, the peak amplitudes 
(i.e., one-half of the peak-to-peak value) required range from 2 to 40 degrees. These can be 
achieved by launching the models in different sabots with different in-sabot tilts. (See, for 
example, Fig. 41.) This requires a number (3 or 4) of different sabot designs. At the larger sabot 
tilts, poor launches can be obtained, because, with the asymmetric sabot structure and 
asymmetric sabot separation, the sabot fingers can re-contact the model after separation, 
causing large, unwanted variations in the pitch and yaw oscillation amplitudes. 

A second way to achieve a range of pitch (or yaw) oscillation amplitudes is to use a single 
sabot design (that orients the model at zero pitch and yaw angles) and “slap” the model with a 
piece of paper or cardstock after the model passes through the sabot stripper. In the NASA 
Ames HFFAF ballistic range, the location of the paper or cardstock is shown in Fig. 23. The 
paper/cardstock is arranged to impact one-half of the model to provide the maximum kick for a 
given paper mass loading (g/cm*). (We use “slap” to describe the action of the paper or 
cardstock impacting the model and “kick” to describe the amount of angular velocity imparted to 
the model by the “slap”.) Figure 42 shows 4 frames from a high speed video showing a LDSD 
model penetrating the slap paper. 
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The paper mass loading (m) to produce a required rms pitch oscillation amplitude can be 
calculated as follows. The linear impulse, I,_, to the model on impacting the paper is given by 


ae = Di mu (42) 


We assume that the impact is a flat, normal impact — no attempt is made to refine the analysis 
by allowing for the paper sliding off the inclined front surface of the model, while maintaining a 
fraction of its initial forward velocity (relative to the model). Since the distance of the center of 
gravity of the semicircle from the axis is 2D/(3z), the angular impulse, I,, to the model on impact 
is 


Lae 
I, =—D*mu 43 
A= Ts (43) 


6 x = (44) 


The angular motion of the model subsequent to impact is given by 


pies. sin( 22") (45) 
P 
where 
0 = instantaneous pitch angle 
Omax = Maximum pitch angle 
t= time 
P = period of pitch oscillations 
Differentiating Eq. (45), we get 
2 2 
0=6... 7 cos( 221 (46) 
P i 
Evaluating Eq. (46) at t = 0, we get 
Z 
Pisa Oi (2) (47) 
P 


Figure 43 shows the rms pitch oscillation angles versus slap paper (or cardstock) mass 
loading for LDSD models with the inflatable stowed (Fig. 11) and deployed (Fig. 6). Two sets of 
data points are shown for the LDSD model with the inflatable deployed, red dots for shots with a 
single slap sheet and blue dots for multiple slap sheets. The purple dots are the data points for 
the LSDS model with the inflatable stowed. It is seen that the theoretical lines agree roughly 
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with the trend liens of the experimental data. For the data with the inflatable deployed, the 
agreement is quite good. For the data with the inflatable stowed, the theoretical line lies up to 
30% above the trend line for the experimental data for the largest angles, but the agreement is 
much better at the smaller angles. 

It has been observed at NASA Ames that stacking sheets of paper or cardstock one 
immediately behind the next can lead to wide variations in the oscillation amplitudes obtained, 
while very consistent results are obtained if single sheets of paper or cardstock of varying 
thicknesses are used. It is believed that these variations are caused by the ejecta from the first 
sheets blowing away part of the succeeding sheets before they can produce their full impact 
momentum. This effect is believed to be very sensitive to the exact arrangement and spacing of 
the sheets. The effect can be seen in Fig. 43, where four of the blue data points, with multiple 
slap sheets, are significantly below the trend line, while the red data points, with single slap 
sheets, group much more tightly around the trend line. All of the purple data points, with single 
slap sheets, group about as well as the red data points. At Ames, the paper/cardstock sheets 
are stacked together without separation. On the other hand, at the Eglin ARF range, multiple 
sheets of paper or film with spacings of several centimeters have been successfully used. 

Figure 44 shows data corresponding to the data of Fig. 43 for a proprietary model launched in 
the NASA Ames HFFAF facility. Figure 44 shows the difference in the performance of the slap 
paper for the supersonic (M = 1.3 — 1.7) and subsonic (M = 0.6) cases. This is likely due, at 


least in part, to the difference in Cn, between the supersonic and subsonic cases. It also shows 


that, for the supersonic cases, when the cutout in the slap paper is used, about 2.2 times the 
paper mass loading is required to achieve the same rms pitch oscillation amplitudes. The semi- 
circular cutout was 1” in diameter and the model diameter was 1.5”. The cutout was used 
because the center of the front of the model was a very thin shell and there was concern that 
the slap paper impact could permanently deform the shell. 

For models with relatively thin, fragile aerobrakes surrounding a compact aeroshell, estimates 
of the stresses produced by the impact of typical slap paper often greatly exceed the strengths 
of the aerobrake materials. One possibility to avoid this problem is to “slap” the model by 
passing it through, say, 30 very light (~O.0006 cm thick) Mylar diaphragms spaced, say, 1 cm 
apart to avoid the ejecta issue discussed above. A second possibility, which provides an even 
smoother “slap” would be to pass the model through a gas jet which is heavier or lighter than 
the surrounding air. The heavier jets (e.g. Ar, COz) would increase the aerodynamic stress on 
the aerobrake, but a lighter (e.g., He) jet would not increase the aerodynamic stresses, but 
would provide a “negative” slap since the aerodynamic forces are unbalanced in the opposite 
direction. Air on the side of the model opposite the helium would be doing the “slapping”. It is 
clear that for both these options, the slapping medium, Mylar sheets or a gas jet, must be 
arranged to cover one half of the model as it flies through the medium. 


Appendix A 


In this appendix, we present the key parameters for the flight trajectory and the ballistic range 
simulation of the LDSD model with the inflatable deployed discussed in Sec. 5 of the main text. 
The list below gives key parameters for the flight vehicle and the normalized radius of gyration 
and the density for the preliminary range model. 


D, flight vehicle diameter = 6 m 

m, flight vehicle mass = 1369.9 kg 

Imin, flight vehicle minimum moment of inertia = 1513.8 kg-m* 
Imax, flight vehicle maximum moment of inertia = 2077.4 kg-m? 
| in the analyses herein is taken to be Imax 

V, flight vehicle volume = 23.017 m° 

Os, flight vehicle mean density = m/V = 59.52 kg/m° 
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JU 
A, flight vehicle frontal area = qo = 28.27 m* 


Xcg/D, normalized axial location of center of gravity, distance aft of vehicle nose = 
0.1583 


I 
fg, flight vehicle radius of gyration “fe = 1.231 m 
—— for flight vehicle = 0.04212 


—— for preliminary range model = 0.05993 


os for preliminary range model = 7800 kg/m” 


Table A1 below gives the trajectory data, the CFD aerodynamics data for the flight vehicle, 
the density ratios p/p, for the flight vehicle and the preliminary model design, the normalized 
pitch oscillation wavelength (A/D) and the K parameter [see Eq. (4) of Sec. 5] for the flight 
vehicle and the preliminary model design. 


Table A1. Trajectory, aerodynamic and other key data for the LDSD vehicle with the 
inflatable deployed and the corresponding range model. 


point 1 point 2 point 3 point 4 point 5 
Po MUTT 8957 | 85O | 8.00] 2479 | 2.365 
Pp | kg/m* | 1.107 x 10" | 1.058 x 10 | 1.088 x 10 | 1.379 x 10" | 1.524 x 10° 
pC Re | | 5.126 x 10° | 4.337 x 10° | 3.821 x 10° | 4.001 x 10° | 4.220 x 10° 
C 


PG 0.1935 0.185 0.175 0.164 0.161 
P 6 Syt~—“<‘i;*;‘C ‘;SCOC‘#*‘«A‘SS4 1.468 1.492 1.530 1.544 


C ] 1.127 1.155 1.192 1.249 1.268 


Ly 
| pep (flight) | | 53,760 | 56,250 | 54,700 | 43,160 | 39,050 
MD | | 853.8 | 870.1 | 875.3 | 344.3 | 330.8 
pom 8141735 | 1710 | 1865 | 1.941 
| K(flight) | | 1.525 | 1560 | 1.604 | 1.657 | 1.673 
_pslp(pre.model) | | 37,730 | 39,480 | 38,390 | 30,290 | 27,410 
| K(pre.model) | | 1.820 | 1862 [| 1914 | 1.977 | 1.997 ___ 


“n is the number of pitch oscillation cycles in the range with a model diameter of 3.56 cm and a 
range length of 22.86 m 
“pre. model” denotes preliminary model 


We note that the number of cycles in the range shown in Table A1 is somewhat less than the 
goal of 2.0 cycles. To alleviate this problem, the pressure in the range was raised somewhat 
above that implied by the p,/p (preliminary model) values. This brings the number of cycles back 
towards the goal of 2.0, while accepting some mismatch (up to 20 - 30%) of the Re and A/D 
values. In the actual test, the number of cycles in the range varied from 1.84 to 2.4. Only 7 of 37 
shots had less than two cycles in the range and only 1 shot had less than 1.9 cycles in the 
range. 
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Appendix B 


In this appendix, for the LDSD vehicle with the inflatable deployed and the range simulation of 
same, we give the key parameters necessary for the calculation of the velocity and Mach 


number drops a and dMrange and the velocity Ax (see Sec. 7). The list below gives some 
U |, 
pertinent parameters for the flight vehicle and the preliminary range model. 


—— for flight vehicle = 0.04212 


—— for preliminary range model = 0.05993 


L, range length = 22.86 m 
D for preliminary range model = 3.56 cm 


Table B1 below gives the data for the calculations for the five data points along the trajectory 
of the LDSD vehicle with the inflatable deployed and simulation of same. 


Table B1. Data for calculation of velocity and Mach number drops and slow down Axyange 
for LDSD vehicle with the inflatable deployed and simulation in the range. 


Variable Units Trajectory Trajectory Trajectory Trajectory Trajectory 

point 1 point 2 point 3 point 4 point 5 

M(average) | | __3.957 2.479 2.365 
C 


¥ i 0.1935 0.185 0.175 0.164 0.161 
6 Spt~<“<‘iU;‘“‘C‘édSCO!;*C‘~ 4S“ 1.468 1.492 1.530 1.544 


| MD | | 353.8 370.1 375.3 344.3 330.8 


du 
Ewe 0.0353 0.0357 0.0378 0.0451 0.0482 
U A, flight 
Hp 0.0444 0.0448 0.0475 0.0566 0.0606 
U A,range 


ee RO LT: 0.272 0.244 0.262 0.278 
120.8 134.5 


Appendix C 


In this appendix, for the LDSD vehicle with the inflatable deployed and the range simulation of 
same, we give the key parameters necessary for the calculation of the maximum swerve in the 
range due to the oscillating lift due to pitch oscillations, ymax [see eq. (26) of Sec. 8]. The list 
below gives some pertinent parameters. 


2 
—~ for preliminary range model = 0.05993 
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D for preliminary range model = 3.56 cm 
Omax, Maximum pitch angle = 10 degrees 


Table C1 below gives the data for the calculations for the five data points along the trajectory 
of the LDSD vehicle with the inflatable deployed and simulation of same. 


Table C1. Data for calculation of maximum swerve due to oscillating lift vector for LDSD 
range model with the inflatable deployed. 


Variable Units Trajectory Trajectory Trajectory Trajectory Trajectory 
point 1 point 2 point 3 point 4 point 5 


—  M | 3.957 2.479 2.365 


ae a 0.1935 0.185 0.175 0.164 0.161 
PG 1.127 1.155 1.192 1.249 1.268 


Ly 
0.253 0.245 0.219 0.201 0.187 


Appendix D 


In this appendix, for the LDSD vehicle with the inflatable deployed and the range simulation of 
same, we give the key parameters necessary for the calculation of the maximum swerve in the 
range due to steady-state lift, y [see Eq. (32) of Sec. 9]. Note that the actual model has zero 
mean lift and hence the actual y values are all zero. A lift coefficient of 0.23 was assumed in this 
appendix to illustrate the calculation technique. The list below gives some pertinent parameters. 


2 
r 
aE for preliminary range model = 0.05993 


D for preliminary range model = 3.56 cm 
L, full length of range from gun muzzle to last shadowgraph station = 32.92 m. 


Table D1 below gives the data for the calculations for the five data points along the trajectory 
of the LDSD vehicle with the inflatable deployed and simulation of same. 


Table D1. Data for calculation of swerve due to (assumed) steady-state lift vector for 
LDSD range model with the inflatable deployed. 


Variable Units Trajectory Trajectory Trajectory Trajectory Trajectory 
point 1 point 2 point 3 point 4 point 5 


3.957 2.479 2.365 


ae 
es ae 0.1935 0.185 0.175 0.164 0.161 
ee 


Since the estimated maximum lift height while still within the viewing area of the shadowgraph 
stations is 47 cm (see Sec. 9, main text), at the lowest Mach number point, we would be just 
about at the limit with C_ = 0.23. For the higher Mach number points, somewhat higher C_ 
values, up to ~0.31 would be acceptable. 
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Fig . 1. Arrangement of gun, dump tank, sabot stripper and shadowgraph range for the NASA Ames HFFAF facility. 
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Fig. 2. Model inside 3 of 4 serrated sabot fingers. Fourth sabot finger is off to the right. 
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Fig. 3. Sketch (not to scale) of representative two-stage light gas gun. 
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Fig. 4. Representative plot of pitch angle of attack versus shadowgraph station number. 
Note lines joining extrema. 
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Fig. 5. LDSD vehicle with inflatable deployed. Fig. 6. Ballistic range model of LDSD vehicle with inflatable deployed. 
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Reynolds number 


Ames HFFAF LDSD simulation 
trajectory with inflatable 
deployed (solid red lines) 

Numbers beside data points 
and lines are Mach numbers. 


Solid blue lines show parameter spaces) 


that can be covered by the Ames 
HFFAF range for M = 3.96 

and max model diameter of 3.56 cm 
and minimum model diameter of 
1.27 cm. 

Dashed blue lines show extension of 
parameter space for max mode! 
diameter of 4.88 cm. 

Dotted biue lines show extensions 
of p, high limit lines 

Dashed red lines show various limits 
which are discussed in the text. 
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Fig. 7. Ames HFFAF LDSD simulation trajectory with inflatable deployed. Solid blue lines show facility simulation capability for Mach 3.96 and a 
maximum model diameter of 3.56 cm. Dashed blue lines show facility simulation capability for a maximum model diameter of 4.88 cm. 
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Ames HFFAF LDSD simulation 
trajectory with inflatable deployed. 
Numbers beside data points 
and lines are Mach numbers. —- LDSD deployed flight, variable models 
Glue Hines show -2- LDSD deployed final (sing! 
parameter spaces that ployed range, final (single) 
can be covered by the Ames model 
HFFAF range for Mach 3.96 
and 2.37 with model diameters 
ranging from 1.27 to 3.56 cm. 
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Density ratio, R = p/p 


Fig. 8. Ames HFFAF LDSD simulation trajectory with inflatable deployed. Red data points show trajectory with variable models. Green data points 
show trajectory with a single model. Solid and dashed blue lines show facility simulation capability for Mach 3.96 and Mach 2.37 and a maximum 


model diameter of 3.56 cm. 
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Ames HFFAF LDSD simulation 
trajectory with inflatable deployed. 
Numbers beside data points 
and lines are Mach numbers. 
Blue lines show parameter spaces 
that can be covered by the Ames 
HFFAF range for Mach 3.96 
and 2.37 with model diameters 
ranging from 1.27 to 3.56 cm. 
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- - LDSD HFFAF deployed limits 
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Fig. 9. Ames HFFAF LDSD simulation trajectory with inflatable deployed. Red data points show trajectory with variable models. Green data points 
show trajectory with a single model. Blue lines show facility simulation capability for Mach 3.96 and Mach 2.37 and a maximum model diameter of 
3.56 cm. Purple line shows points at Mach 3.0 for a single model obtained by varying the range pressure. 
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Fig. 11. Ballistic range model of LDSD vehicle with inflatable stowed. 


Fig. 10. LDSD vehicle with inflatable stowed. 
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Ames HFFAF LDSD simulation 
trajectory with inflatable stowed. 
Numbers beside data points 


and lines are Mach numbers. 
Blue lines show parameter spaces that —+- LDSD stowed flight, variable models 


can be covered by the Ames += LDSD stowed range, final model 
HFFAF range for Mach 3.45 

and model diameters ranging from 

1.27 to 3.56 cm. 


LAMPCALO5 
Graph 1242 


Reynolds number 
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Fig. 12. Ames HFFAF LDSD simulation trajectory with inflatable stowed. Red data point is un-attainable (in range) condition. Upper green data point 
is closest attainable range condition. Solid blue lines show facility simulation capability for Mach 3.45 and a maximum model diameter of 3.56 cm. 


Lower green data point is alternate range operating condition at Mach 3.15. 
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Fig. 13. HEART vehicle. 
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Ames HFFAF HEART simulation trajectory 

Numbers beside data points are Mach numbers. 

Assumed minimum of 6 data points/per cycle, max 
mode! diameter equals 0.80 times barrel diameter 

Blue and green lines show parameter spaces that can 
be covered by the Ames HFFAF, for M = 2.0 and 


M = 0.8, respectively For 61 mm 
gun 


5 6789 : 
10 
Density ratio, R = p/p 


Fig. 14. Ames HFFAF HEART simulation trajectory. Blue lines show facility simulation capability for Mach 2.0 and maximum model diameters of 
4.88 and 3.56 cm. (61 mm and 44 mm guns). Green lines show corresponding capabilities for Mach 0.8. 
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Fig. 15. Capabilities of the Ames GDF range for the HEART vehicle at the Mach 0.8 flight data point, using the Ames 44 mm powder gun. 
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Ames HFFAF HEART simulation trajectory 

Numbers beside data points and lines are Mach numbers. 

Assumed minimum of 6 data points/per cycle, max 
mode! diameter equals 0.80 times barrel diameter 

Blue lines show parameter spaces that can be covered 
by the Ames HFFAF, for M = 2.0 

Green lines show parameter spaces that can be covered 
by the Ames GDF, for M = 2.0 

Purple lines and square show parameter spaces that 
can be covered by the DRDC range, for M = 2.0. 

Lines for 155 mm gun; square for 110 
mm gun. 
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Fig. 16. HEART simulation trajectory. Blue lines show facility simulation capability for Ames HFFAF facility for Mach 2.0 and maximum model 
diameters of 4.88 and 3.56 cm. (61 mm and 44 mm guns). Green lines show the corresponding simulation capability for the Ames GDF facility. 
Purple lines show the corresponding simulation capability for the DRDC facility with maximum model diameters of 8.80 and 12.4 cm (110 mm and 


155 mm guns). 
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Ames HFFAF HEART simulation trajectory 


Numbers beside data points and lines are Mach numbers. 


Assumed minimum of 6 data points/per cycle, max 
mode! diameter equals 0.80 times barrel diameter 

Blue lines show parameter spaces that can be covered 
by the Ames HFFAF, for M = 0.8 

Green lines show parameter spaces that can be covered 
by the Ames GDF, for M = 0.8 

Purple lines and square show parameter spaces that 
can be covered by the DRDC range, for M = 0.8. 
Lines for 155 mm gun; square for 110 mm gun. 
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Fig. 17. HEART simulation trajectory. Blue lines show facility simulation capability for Ames HFFAF facility for Mach 0.8 and maximum model 
diameters of 4.88 and 3.56 cm. (61 mm and 44 mm guns). Green lines show the corresponding simulation capability for the Ames GDF facility. 
Purple lines show the corresponding simulation capability for the DRDC facility with maximum model diameters of 8.80 and 12.4 cm (110 mm and 


155 mm guns). 
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SIDE VIEW (SECTION) 
Fig. 19. Ballistic range model with sideways press-fit weight and “anti-weight” to adjust lateral position of CG. 
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Fig. 21. Three piece tabbed model. Several methods of adjusting the axial and lateral position of the CG of this model are discussed in the text. 
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Fig. 22. Model inside 3 of 4 serrated sabot fingers. Fourth sabot finger is off to the right. 
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Fig. 23. Arrangement of gun, dump tank, sabot stripper and shadowgraph range for the NASA Ames HFFAF facility. 
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Fig. 24. Launch package designed to use both aerodynamic (ram air) and muzzle blast sabot separation forces. 
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(a) x= 1.086 m (b) x =1.239 m 


(c) x =1.990 m (d) x=2.167m 


Fig. 25. Frames from sabot separation videos for the launch package shown in Fig. 21. Distances are from the muzzle of the gun. 
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(a) x=1.069m (b) x=1.215m 


(c) x=1.979 m (d) x= 2.144m 


Fig. 26. Frames from sabot separation videos for a launch package identical to that shown in Fig. 21 except without the muzzle blast cavity. 
Distances are from the muzzle of the gun. 
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MODEL 


Fig. 28. Launch package with obturator cup (beneath sabot fingers) designed for aerodynamic sabot finger separation. 
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Fig, 37. X-ray pictures showing successful launch (upper picture) and failed launch (lower picture) of 5 degree cones. Failed launch is at a higher 
velocity and with a thinner pusher plate. In the failed launch, the model “fell” through the pusher plate and “cookie cut” out the pusher plate cone. 
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Fig. 38. Aerodynamic forces on free flying model with compact centerbody surrounded by relatively thin aerobrake. (a) is one piece model and (b) Is 
two piece model 
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Fig, 42. Four frames from a high speed video showing a LDSD model penetrating the slap paper. 
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Fig. 43. RMS pitch oscillation amplitudes versus slap paper mass loading for two LDSD models 


launched in the NASA Ames HFFAF facility. 
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Fig. 44. RMS pitch oscillation amplitudes versus slap paper mass loading for a model 
launched in the NASA Ames HFFAF facility. 
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